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ABSTRACT 

We present a study of the Ha line emission from a sample of 1482 optically- 
selected, morphologically-classified bright galaxies (median redshift of 0.05) de- 
rived from the Sloan Digital Sky Survey. The luminosity function is calculated 
for each morphological class and for the total sample. The luminosity function 
fitted with the Schechter form gives a slope a = — 1.43±0.10 for the total sample 
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and the Ha luminosity density is 10 ' ' -o orh erg s~ Mpc~ , where the first 
error is statistical and the second is systematic. This value is consistent with 
that derived by Gallego et al. (1995), but this agreement is caused by a fortu- 
itous cancellation of their neglect of stellar absorption that affects the estimate 
of extinction corrections and a significant sample incompleteness of emission line 
galaxies. The fraction of Ha emitters monotonically increases from early (a few % 
for ellipticals) to late types (100% for irregular galaxies), whereas strong emitters 
exist in all classes of morphological types. We find that 83% of the luminosity 
density comes from spiral galaxies, 5% from irregular galaxies, and 9% from early 
type galaxies; a small number of morphologically disturbed galaxies contribute 
by 3%. 

Subject headings: galaxies: fundamental parameters 



1. Introduction 

Much attention has been devoted to the global star formation rate of galaxies as a 
function of redshift. These studies have advanced our understanding of how galaxies evolved 
towards the present epoch (Gallego et al. 1995; Madau et al. 1996; Lilly et al. 1996; 
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Treyer et al. 1998; Tresse & Maddox 1998; Steidel et al. 1999; Glazebrook et al. 1999; 
Sullivan et al. 2000; Tresse et al. 2002). Despite its obvious importance, the estimate of 
the star formation rate at zero redshift, which is the boundary value for the star formation 
history, still largely relics on the pioneering work of Gallcgo ct al. (1995) which was based 
on Schmidt objective prism plates searching for strong Ha emission galaxies with observed 
equivalent widths EW(Hq;+[NII])>10A. The evolution of galaxies is inferred by comparing 
the global star formation rates at higher redshift with Gallego et al.'s value a,t z ^ 0. Radio 
investigations have suggested star formation rates higher than Gallego et al.'s by a factor 
of two (Cram 1998; Serjeant, Gruppioni & Oliver 2002). The follow-up work for Gallego et 
al.'s sample by Gil de Paz et al. (2000) gave valuable information as to the nature of star 
forming galaxies, yet the question of when galaxies show star formation activity among their 
entire population remains uncovered. Also a poorly understood aspect is the star formation 
activity as a function of morphology of galaxies, which constitutes another dimension of 
galaxy formation. 

In fact, there are a few pieces of work which studied the star formation rate as a function 
of morphologies of galaxies (Kennicutt & Kent 1983; Ryder & Dopita 1994; hereinafter 
RD94). These studies, however, were based on nearby galaxy samples, which were preselected 
by choice, while ideally one want to study homogeneous galaxy samples. The primary 
difficulty of this approach for nearby galaxies is to obtain accurate spectroscopic information 
for a homogeneous galaxy sample over a large area of the sky. The Sloan Digital Sky 
Survey (SDSS; York et al. 2000) overcome this difficulty with its spectroscopic data for a 
large number of homogeneously selected galaxies. The SDSS data also allow significantly 
more accurate treatments as to the control over selection effects, such as those expected in 
line-selected samples, and systematic errors arising from stellar absorption, extinction, and 
contamination from AGN. 

In this paper we use the optically-selected and morphology-classified bright galaxy sam- 
ple from the Sloan Digital Sky Survey for the northern equatorial stripe. The sample is 
flux-limited with r* < 15.9 mag after dcrcddening corrections, and morphological classifica- 
tion is given by visual inspections. This sample, with median redshift z = 0.05, has been 
used to derive luminosity functions (LF) of morphologically selected galaxies (Nakamura et 
al. 2003; hereinafter N03). The SDSS spectroscopy covers from 3900 to 9000 A with the 
resolution of 2-5 A. Redshifts are identified and EWs are measured with the spectroscopic 
pipeline. We use Ho; as an indicator of star formation activity. The SDSS spectroscopy is 
sufficient to resolve Ha and [N II]AA6548,6583 doublet. The advantage of using Ha as the 
star formation indicator is that this emission arises from the ionising flux that is produced 
by massive short-lived stars compared, for instance, to the indicators using the UV light at 
1500-2800A which survives, or even increases, for a significant duration of time after stars 
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formed (e.g., Tresse & Maddox 1998; Glazebrook et al. 1999). Another advantage is that 
the extinction correction is smaller and readily corrected for by using the Balmer line ratio. 
This contrasts to UV emission, which suffers from large extinction but also an uncertainty 
in extinction laws (Calzetti et al. 1994). The most important disadvantage with the SDSS 
spectroscopy for our purpose is that it uses fibres of small apertures, and the aperture cor- 
rection is essential. We circumvent this problem by using the empirical radial profile of Ha 
obtained by RD94 to calculate the total Ha flux. 

Our prime purpose is to investigate the Ha cmissivity distribution for each class of 
morphological type of galaxies, and its contribution to global Ha luminosity in the nearby 
universe {z < 0.12) using a well-defined optically fiux-limited sample. We prefer to use 
traditional morphological classification obtained by visual inspections, since spectroscopic or 
colour indicators are too sensitive to small amounts of star formation activity, whereas visual 
morphology reflects dynamical evolution of galaxies more faithfully. The drawback of this 
approach is the small size of the sample (1600 galaxies). Notwithstanding, we emphasise that 
the present study is based on a galaxy sample with homogeneous morphological classification 
and accurate photometry; the resulting Ha emission sample is significantly larger than those 
used in previous studies. 

In section 2 we describe the data and various corrections that are crucial to our study. 
The Ha LF is derived in section 3 for each morphological class. A detailed comparison is 
given with the work of Gallcgo et al. (1995) for the total Ha luminosity density in this 
section. Morphological-typc and colour dependence of Ha emissivity is presented in section 
4, and implications of our results are discussed in section 5. 

We adopt the cosmology Q = 0.3 and A = 0.7 throughout this paper, although the 
physical parameters of our objects, being at such low redshift, depend little on Qq and Aq. 
We use the standard notation h = ifo/lOOkm s~^Mpc~^ to demonstrate the dependence of 
quantities upon Hq. 

2. The emission line galaxy sample 
2.1. Data 

Our sample is derived from the photometric images (Gunn et al. 1998; Hogg et al. 
2001; Pier et al. 2003) of Run 752/756 of early SDSS observations [the area lies along the 
celestial equator for 9''40"'36" < a(J2000) < 15''43"*53^]. The observations cover 229.7 
deg^ in the five colour band photometric system (Pukugita et al. 1996), calibrated by the 
standard star network established by the work at USNO (Smith et al. 2002). The bulk of 
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data were published in Early Data Release (Stoughton et al. 2002), and are also included in 
Data Release One (DRl, Abazajian et al. 2003). Spectroscopic information is obtained for 
galaxies that satisfy accurately defined selection criteria (Strauss et al. 2002; Blanton et al. 
2003b). The morphological classification was performed by visual inspections for all galaxies 
with dereddened Petrosian magnitude of r*p < 15.9 included in the rectangular area of Run 
752/756^. The morphological index T = — 6 is assigned to E, SO, Sa, Sb, Sc, Sd and Im. 
The galaxies that show largely disturbed morphologies are given T — —1. The size of galaxy 
sample with spectroscopic information is 1600 (See Table 1 of N03). 

Among 1600 galaxies, wc drop 21 objects because of cither low (< 85%) rcdshift confi- 
dence or misallocation of fibres to overlapping stars, 29 galaxies because of poor photometry 
due to deblending of neighboring objects, and 59 galaxies that fall either z < 0.01 or 2; > 0.12. 
We remove an additional 9 galaxies with r*p < 13.2, at which magnitude the spectroscopic 
survey completeness is lowered due to the SDSS target selection algorithms. These cuts leave 
1482 galaxies for our consideration. This data set is the same as the sample from which the 
morphologically-dependent luminosity function was derived in N03 (sample 5 in Table 1 of 
N03), and is recapitulated in Table 1 below. These selections do not introduce any significant 
bias into the sample with respect to morphology or luminosity. 

The SDSS spectroscopy was carried out using fibres of 3" diameter placed on the centre 
of each object with F/5 optics. The 2-5A spectral resolution of the double spectrographs 
resolves Ha emission from [N11]A6548 and [N11]A6583. This contrasts to earlier works where 
Ha is often blended with nitrogen lines (e,g., Sullivan et al. 2000; Trcssc & Maddox 1998). 
Redshifts and EWs are automatically measured by the SDSS spectroscopic pipeline and 
are catalogued. We use redshifts converted into the values in the Galactic standard of rest 
according to de Vaucouleurs et al. (1991). There is httle confusion in the measurement of 
Ha equivalent widths: visual inspection of spectra shows that emission lines with greater 
than Ri0.5A in the observed frame are accurately measured for the majority of galaxies. 
Manual remeasurements of EWs for selected galaxies using IRAF show that the error is less 
than O.5A for small EWs and less than 10% for larger EWs. The error mostly arises from 
the fit to the continuum. We shall work with the rest-frame equivalent width (hereafter we 
always refer to the rest-frame EW unless otherwise noted) with measurements at more than 
2.5(7. A total of 30 galaxies have EW> lA yet the detection confidence is less than 2.5a: 
those galaxies are omitted from our main sample. We find that most cases this occurs when 
the Ha fine is too broad and/or emission of [Nil] is much stronger than that of Ha, and 



Although new photometry is available in DRl, we use old photometry of EDR since the morphologically 
classified catalogue is made using the sample based on old photometry and the completeness would be lost 
in faint magnitudes if new photometry is adopted. 
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automated deblending of the Ha lines is not properly executed. So, those that are missed 
by 2.5(7 criterion for EW> lA are AGN like galaxies. Finally, we dropped one galaxy for 
which the fibre is centred on a HII region off the centre of galaxy. Thus, our emission line 
galaxy sample contains 665 galaxies. 

We show in Figure 1 the EW distribution of Ha for morphologically-classified sample, 
from the latest type (top panel) to earlier types. The bottom panel is the total sample. 
A negative EW signifies absorption. The solid lines correspond to our main sample with 
EW> lA and the detection at > 2.5a, and dotted are galaxies that do not satisfy these 
criteria. Note that the ordinate is a logarithmic. The mean of emission line strengths for 
each morphological class is indicated by the diamonds. In this average (indicated by dotted 
histogram), we include galaxies with EW< lA and with lower confidence detection, as well 
as those that do not show emission at all. AGN-like galaxies, as discussed in Sect. 2.4, are 
denoted with shading. The numbers are presented in Table 1, where emission line galaxies 
are listed separately for non-AGN and AGN galaxies. The bottom panel shows that for 
EW> SA the galaxies that are excluded from the sample by 2.5a detection are only a few, 
and the sample is virtually complete. The sample incompleteness for EW= 1 — sA is about 
14%. We also try to include those galaxies which do not satisfy the > lA or the 2.5a 
criterion in our analysis and examine how the results change. 



2.2. Stellar Absorption Corrections 

The measured Balmer emission line equivalent widths must be corrected for stellar 
absorption. The absorption not only decreases Ha emission fiux, but also significantly affects 
estimates of reddening corrections inferred from Balmer line ratios. Figure 2 shows the 
relation between EW(Ha) and EW(H/9)xFc(H/3)/Fc(Ha), where Fc is the continuum fiux 
densities at the two lines [Fc(H/?)/Fc(Ha) ~ 1.0 for galaxies with EW(Ha)< 25A in our 
sample]. The relation plotted in this figure is equivalent to that between the line fiuxes, 
F(Ha) and F(H/3). Considering galaxies with EW(Ha)>10.0A, for which corresponding 
EW(H/3) is mostly larger than I.OA and hence the H/3 measurement is reliable, we obtain 
with a fit, 

EW(H/3) X ^'^^fl = -1.58 + 0.246 x EW(Ha), (1) 

Fc(Ha) 

where 3a rejection is applied. The existence of a non-zero constant term indicates the 
presence of stellar absorption. If we take the Ha absorption to have an equivalent width 
of — l.sA, as indicated from the conspicuous clustering of non-emission line galaxies seen in 
the figure, we obtain from eq. (1) EW(H/3)abs = — 2.0A, which is also consistent with the 
position of EW(H/3) for non-emission fine galaxies. 
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An independent check of the effect of stellar absorption is carried out assuming EW(H7)abs 
=EW(H/5)abs (e.g., McCall, Rybski & Shields 1985) for 168 galaxies with EW(H7) de- 
tection at > 2.5(T out of 665 emission line galaxies. After extinction corrections and as- 
suming the standard Balmer line ratios, as we describe in the next subsection, we obtain 
EW(H7)abs ~ — I.5A. (see Figure 3) independent of EW(Hq;). 

In our following analysis wc take EWabs(HQ;)=-1.8A and EWabs(H/9)=-2.0A for all 
galaxies. We note that the estimate of extinction corrections depends little on the absolute 
value of the stellar absorption in so far as we use the H/3/Hq; line ratio from (1). 



2.3. Extinction Corrections 

The extinction correction for the line flux is estimated from the Hq;/H/3 Balmer line 
ratio, which is written 

^i^*^) _ r, 1 n-c[/(Ha)-/(H/3)] (o) 

where D is the ratio of the intrinsic fluxes emitted in the nebula, /(A) = A;(A)/A';(H/?), and c — 
OAk{}ip)E{B - V) with k{\) the Whitford extinction curve that satisfies A;(A) = Ax/E{B - 
V). We take Seaton's (1979) law for A;(A), which gives A;(H/5) = 3.68 and k{Y{(5) / k{Y{a) = 
1.48. We also evaluate the extinction using O'Donnell (1994) extinction law with it! = 
3.1 ± 0.2 to investigate the systematic error. We set D — 2.86 assuming Baker- Menzel case 
B recombination for the electron temperature — lO^K and the density Ug — 100 cm~^ 
(Brocklehurst 1971). The change of the temperature by a factor of two modifies D only by 
5%; the change of Ue affects D very httle. 

For galaxies with > 2.5a detection of H/3 [EW(Hq;)>10.0A], eq. (2) is applied to each 
galaxy to derive Ahq. The distribution of extinction is shown in Figure 4 after converting 
it into V band extinction Ay using Av/Affa — 1.28. The mean value of extinction is 
(Ay) = I.IO. This distribution is consistent with that obtained by Sullivan et al. (2000) 
for (z) = 0.2 galaxies {{Ay) = 0.97). The mean extinction also agrees with the value 
derived by Kennicutt (1983). If stellar absorption is not taken into account, the extinction 
correction would become erroneously large by 70%, as noted by Sullivan et al. (2000), and 
the resulting extinction correction would also show an anticorrelation with the Ha equivalent 
width. For rather weak emitters, i.e., EW(Hq;)<10.0A or the H/3 detection below 2.5(7, we 
adopt Ajia—O.Sl {{Ay) — 1.03) which is determined by the Hq;/H/3 fine ratio calculated from 
eq. (1). 
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2.4. AGNs 



Since we are primarily interested in the star formation activity, we must remove Ha con- 
tributions from AGN. We use the hne diagnostic diagram employing [NII]A6583/Hq; versus 
[OIII]A5007/H/3 to identify AGN (Baldwin, Philhps, & Terlevich 1981; Veilleux & Oster- 
brock 1987; see also Kewley et al. 2001). We require that four relevant hnes be detected 
at more than 2.5a. Then, 363 galaxies out of 665 galaxies pass this criterion. The line 
diagnostic we use is immune to dust extinction, since the lines of each pair are close enough 
in wavelength. 

We use the criterion set by Kauffmann et al. (2003), which was apphed to the SDSS 
galaxy sample to select AGN: 



, JOIII1A5007, 0.61 

where Ha and H/3 should be corrected for stellar absorption. 89 galaxies satisfy this criterion 
and are identified as Seyfert II or LINER. The frequency of AGN thus identified is 6%, 

compared to 18% in the Kauffmann et al. sample. The difference arises from weak H/3 
line emitters, for which special effort is made by Kauffmann et al. to measure weak lines 
with higher signal-to-noise ratio. In our sample about 3% of galaxies with B.f3 emission of 
EW>1A are missed by our 2.5a criterion, but this fraction increases to 26% for EW=0.5- 
lA and rapidly increases for EW< O.SA . Most of the difference in numbers of AGN between 
Kauffmann et al. and this study arises from galaxies with EW(H/3) < 0.5A, which are not 
important to us. If we use two-hne diagnostic with only Ha and [Nil] [log([NII] A6583/Ha) > 
—0.3], the number of AGN increases to 235 galaxies which is 16% of the total number of 
galaxies. In the next section we derive the Ha luminosity functions, the result of which only 
weakly depends on whether we reject 89 AGN or 235 'AGN'. 



2.5. Aperture corrections 

Since the line flux is measured with fibres of a 3" aperture, aperture corrections are 
essential to estimate total Ha luminosity of galaxies. There is also a problem with spec- 
trophotometric calibration, the error of which is about 20% (Stoughton et al. 2002), 
which we confirmed for our sample. 

The aperture correction is made by using empirical radial profiles of Ha and broad band 
fluxes obtained by RD94. These authors studied Ha and / band surface brightnesses for 34 
nearby SO and spiral galaxies. We show in Figure 5 the Ha flux versus i' band flux (7 mag 
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of RD94 is translated to i! system using Pukugita et al. 1995). The dashed (dotted) hnes 
show the trace of Ha and i' band fluxes of RD94 photometry: it moves left /upwards as the 
aperture increases and the circle denoting their 'total' flux. Superimposed in the bundle of 
traces are SDSS spcctrophotomctric data of Ha and i*-synthetic flux in the 3" aperture. We 
sec that the SDSS data overlap very well with the RD94 traces, except for several galaxies 
(denoted by dots) which deviates from the SDSS data points. 

We construct a composite growth curve of the ratio of Ha to / band flux from 28 
galaxies, rejecting 6 galaxies that are outliers of the SDSS data. We integrate the Ha and / 
band fluxes from the centre of galaxies outwards, and the integrated Ha fluxes are evaluated 
against the integrated / band fluxes, with both total fluxes normalised to unity. We do not 
observe systematic trends in growth curves that depend on morphologies of galaxies^. To 
obtain the broad band Petrosian flux defined by the SDSS, we extrapolate the measured 
flux assuming an exponential profile with the scale length determined by RD94 (this gives 
~ 8 — 20% correction). For the Ha flux, there is a natural cutoff at //Ha ~ 2.5 x 10^^ erg 
s~^kpc~^, as shown by Kennicutt (1998), and the limit of RD94 measurement agrees with 
this cutoff; the integration of the flux measured by RD94 gives the total Ha flux. The growth 
curve thus constructed is shown in Figure 6, together with those of individual galaxies (the 
error bars show one sigma at representative points). The composite growth curve shows that 
Ha emission activity is somewhat more active in the outer region of galaxies than in the 
central part, although for individual galaxies the growth curves show a variety. Note that 
our 3" aperture corresponds to (2.1 — 3.9)/i~^ kpc z — 0.05 — 0.10. 

We flrst renormalisc the spectrophotometric flux by adjusting i' band synthetic magni- 
tude obtained from spectrophotometry to the corresponding 3" aperture photometric magni- 
tude to remove the calibration error of spectrophotometry, and then apply the growth curve 
to find the total Ha fiux from its spectrophotometric fiux and i' band photometric fiuxes at 
the 3" and the Petrosian apertures. Our calculation takes into account the aspect ratio of 
the galaxy images, which gives only a small correction. 

The aperture correction amounts to 2.40±0.83 mag for z ~ 0.015, 2.14±0.59 mag at z ~ 
0.05 and 1.56 ±0.44 mag at z ~ 0.1 (the error stands for rms). The total Ha flux we estimate 
is shown in Figure 7. We do not see any systematic variations of the high luminosity edge 
as a function of redshift, suggesting that aperture corrections are properly done. Systematic 
errors from this aperture correction are estimated by changing the composite growth curve 



'''This docs not mean that the ratio of the total Ha to the total / band fluxes does not depend on 
morphology. RD94 showed that such a ratio actually depends on morphology, whereas the growth curve is 
nearly universal. 
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by ±1(7. If we would apply aperture correction simply by scaling the synthetic Ha flux with 
the broad band fluxes measured at two apertures, we underestimate the total Ha flux by 
20-25%. 

3. The Ha luminosity function 
3.1. Calculation 

We calculate the Ha luminosity function (LF) for five morphological classes and for the 
total sample. Wc treat E and E/SO-SO separately, but when we compute LF we add the 
two together since the number of Ha emission E galaxies (Fukugita et al. 2003) is too small 
to derive a reliable LF for this type. The number of galaxies in each morphological type 
is given in Table 1 above. We reject 89 AGN from the sample. We employ the step-wise 
maximum likehhood estimator of Mobasher, Sharpies, and Ellis (1993), using the r*p band 
LF calculated by the maximum hkelihood estimator in N03. The step-wise Ha LF is given 

by 

,(Lh.)ALh.^ E (4) 

where 0(M**^) is the rp band LF, n{M^*^) is the number of objects that emit Ha luminosity 
between {Lna — AL/2) and {Lna + AL/2) and whose r' band luminosity is in the M^, 
magnitude bin. The function in the sum corresponds to the inverse of a survey vohime 
in a given magnitude bin given by the maximum likelihood estimator. This estimate is 
not affected by inhomogeneity of the sample, unlike the 1/V{max) estimator for which the 
function in the sum is simply the inverse volume of visibility. The results are given in Figure 
8. 

The Ha LF we obtained shows that the Schechter (1976) function 

(f>{Lna)dLua = 0* exp (5) 

gives a reasonable fit. We carried out a fit, rejecting data in the bins where the number 
of galaxies is one (we include such data, however for the LF of Im type, as this entire sample 
is very small). The Schechter function fit suffers a binning artefact. We, therefore, carry out 
10 fits shifting the binning by 1/10 the interval of bins, and smear the resulting parameters. 
We also removed the data in the faintest bin, which may be affected by the binning. The 
results are presented in Figure 8 and the parameters obtained are shown in Figure 9, where 
the error contours are one and two standard deviations of the fit- The parameters and 
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one standard deviation errors are also given in Table 2. We attempt to estimate errors by 
the jackknife method dividing the sky region into 20 segments. We find that those errors are 
comparable to one standard deviation errors we quoted. 

We calculate Ha luminosity density for each morphological class by integrating the 
Schechter function to zero luminosity. This is compared to the luminosity density directly 
calculated as log£(Ha) = Li/V^^{e^) where ^^'"^'^(eff) = n{M'^,^) / (j){M'^,^) that appears 
in eq. (4). The two estimates show good agreement for late types, but the two differ by 
0.1 dcx for E-SO and Im, for which statistics arc very small and the Schechter function fits 
arc rather poor: for these cases the direct sum is probably more reliable, and the luminosity 
densities from the direct sum arc adopted for our discussion. The total luminosity density 
from the Schechter function is log£(HQ;)(erg s~^Mpc"^) = (39.31 ± 0.04) + log/t, which 
agrees with 39.33+log h from the direct sum. The sum of the component luminosity densities 
amounts to 98% of the total (when we discuss the composition of the luminosity densities 
we renormalise the component sum to 100%). Note that the r' band luminosity density of 
N03 is 2.00 X lO^L0/i(Mpc)-^ 

We note that the incompleteness of spectroscopic and/or photometric samples is already 
corrected by the use of r' band LF. The final correction should account for the under or 
overdensity of the bright galaxies in the northern equatorial stripes. The r' band luminosity 
density of N03 has virtually no offset against the global value of Blanton et al. (2003a)^. 

Before discussing the results we address the problem of systematic uncertainties. The 
first issue is the effect of the threshold set for Ha detection (EW(Hq;) >iA). In order to 
test the effect we decrease the detection threshold to EW(Hq;)> —0.8(A) removing the 2.5a 
detection criterion. The new threshold corresponds to selecting galaxies with EW(HQ;)cm > 
1.0(A) after absorption correction. With this change the numbers of galaxies increase to 51, 
158, 394, 337, and 10 for E, SO, early spiral, late spiral, and Im galaxies, respectively (for 
weak line emitters, we do not reject galaxies that possibly satisfy the AGN criterion). While 
this change increases the number of emission line galaxies by 50%, it does not appreciably 
affect the Ha LF for the late-type and the total samples. The change is observed in the 
fainter bins for E and SO galaxies, which are dominated by non or weak Ha emitters around 
the threshold. The brighter part of the E plus SO luminosity function is unchanged. The 
effect on the Ha luminosity density is only +4%. 



^In N03 wc noted that the luminosity density of bright (r* < 15.9) galaxy sample is lower than the 
total spectroscopic sample with r* < 17.88 (Blanton et al. 2001) by 29%. A further study (Blanton et 
al. 2003a), however, showed that what is overdense is the spectroscopic sample of the northern equatorial 
strip compared to the global value. In the r' pass band, the global luminosity density is approximately 
2.0 X lO«X0/i(Mpc)-3. 
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There is some uncertainty in the estimate of the stellar absorption correction. A 
±0.5A uncertainty in the stellar absorption strength yields a —3 to +7% change in the 
luminosity density. A complete neglect of stellar absorption would affect the B.f3/B.a flux 
ratio for weak emitters, and leads to a significant overestimate of extinction corrections: this 
results in the Ha luminosity about 0.4 dex brighter, and the luminosity density 50% larger 
than the true value (NB: 0* decreases). The change of extinction law from the Seaton (1979) 
to the O'Donnell (1994) relations (with R — 3.1) does not significantly affect the results: 
the luminosity density increases by only +2%. Allowing for it! = 3.1 ± 0.2, the uncertainty 
is —2% to +10%. We also allocate an uncertainty of ±10% from the error for the estimate 
of (Ana)- 

If we exclude additional 146 weak H/3 'AGN' as discussed in sect. 2.4, the luminosity 
density decreases by 5%. The shape of the Ha LF changes very little. A complete inclusion 
of all AGN would increase the Ha luminosity density, but it is no more than by +18%. 

The systematic uncertainties we considered are summarised in Table 3, expressed in the 
form of the change in the luminosity density. We also include the cases which, we believe, 
are unrealistic in order to show the effect if we ignore the relevant considerations, as are 
occasionally seen in some previous works. The systematic error is obtained by adding the 
entries with asterisk in quadrature: our estimate is +0.10 dex and —0.07 dex for C 



3.2. The results 

The characteristics of our Ha LF are as follows. 

1. The Ha LF for the total sample is fitted well by a Schcchtcr function with a = 
— 1.43 ± 0.10 and logL* = 41.68 ± 0.10. The Ha luminosity density derived from our LF, 

log £(Ha) (erg s'^Mpc"^) = 39.31 ± 0.04 + log h. (6) 

We estimate the systematic error to be t.olj- 

2. Early- (1.5 < T < 3) and late-type (3.5 < T < 5) spirals are the dominant contrib- 
utors to the total Ha LF. The Ha LFs for these two types of galaxies are similar, and they 
determine the shape of the total LF. The faint end slopes are somewhat steeper than those 
for the r' band LF. Spiral galaxies (1.5 < T < 5) produce 83% of the total Ha luminosity 
density. 

3. The contribution of early type galaxies (E and SO) is only 9% of the total luminosity 
density. In particular the luminosity density from E (8 galaxies) is only 6% that of all early 
type (E to SO) galaxies. The characteristic luminosity of the Ha LF for the E-SO types 
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is about 0.7 dex fainter than those for spiral galaxies. While our 'best fit' to the Ha LF 
indicates a decline towards the faint end, we are not able to conclude whether the Ha LF 
for early type galaxies actually declines because of the small number of galaxies and the 
omission of very weak Ha emitters. The inclusion of weak emitters [EW(Ha)> —0.8(A)] 
lifts the faint end slope, although a is still larger than —1 (a ~ —0.5). 

4. The Ha LF for Im type galaxies shows a steep faint end slope. This result is similar 
to that found for the r' band LF. Im galaxies contribute 5% of the local Ha luminosity 
density. This is rather a significant fraction if we consider the fact that only 0.7% of galaxies 
(10 galaxies) in our sample arc Im type and their broad band luminosity is fainter than that 
of normal galaxies. We also note in Figure 1 that all Im galaxies arc strong Ha emitters. 

5. Most of the galaxies with unclassified morphology (6 out of 7 galaxies) are also strong 
Ha emitters [EW(Ha) > 20A]. These unclassified galaxies contain mergers. 

3.3. Compcirison with Gallego et al. 

Since the work of Gallego et al. (1995) provides the current fiducial value for the star 
formation rate at zero redshift, we shall present a detailed comparison of their results and 
our Ha LF. 

The luminosity function and luminosity density we derived are consistent with the find- 
ings of Gallego et al. who derived a = —1.3 ± 0.2, and logL* = 41.56 ± 0.08 — 2\ogh. 
While this L* looks somewhat smaller than our value, L* and a are correlated such that 
A(logL) ^ -0.4Aa; so if we adjust a = -1.43, then we have logL* = 41.61 ±0.08 - 2 log /i 
which is smaller than our logL* = 41.68 ± 0.10 — 2 log /i only by 1 a. The luminosity density 
log£(Ha) = 39.38 ± 0.04 + \ogh is larger than (6) by about 20% (0.07 dex). We consider, 
however, that these rather good agreements are fortuitous. 

Gallego et al. did not take account of stellar absorption. The neglect of stellar absorp- 
tions leads to a substantial overestimate of extinction corrections (see section 2.3 above). 
We find that their selective extinction E{B — V) tabulated in the catalogue of Gallego et al. 
(1996) show a conspicuous decrease as the Ha equivalent width increases, which is the fea- 
ture that occurs if stellar absorption is ignored. In fact their mean value of E{B — T/) ~ 0.64 
{Ay — 2.05) is 1.9 times as large as the extinction we obtained. This overestimate of the 
absorption makes the luminosity density larger by 0.16 dex than the true value (sec Table 
3). The inclusion of AGN also leads to an overestimate of the luminosity density by < 0.07 
dex. 
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Gallego et al. set the threshold EW(HQ;+[NII])>10A(effective Ha EW 7A); this selec- 
tion would cause a 10% decrease in the Ha luminosity density. What gives a strong effect is 
their sample incompleteness, which compensates the significant overdensity seen above. The 
surface density of Ha emitters with EW(Ha+[Nll])>10A is 176/471=0.37 per sq. deg. This 
is compared with 306/230=1.33 per sq. deg. when wc select emission line galaxies from the 
SDSS sample in a manner that satisfies the same EW threshold and z < 0.045 as in Gallego 
et al. The EW distributions of the two samples are compared in Figure 10, which shows 
where Gallego et al.'s sample is incomplete. Note that we have applied aperture corrections 
for the EW distribution of the SDSS sample so that the line flux represents that from the 
entire galaxy. It appears that the underestimate due to sample incompleteness cancels the 
overestimates from the neglects of stellar absorptions and contamination from AGN. 



4. Correlation of Ha emissivity with galcLxy properties 

4.1. Ha emissivity and galaxy morphology 

Figure 1 and Table 1 above show the clear trend that the fraction of Ha emission 
galaxies increases as T moves from early to late types. If we set the threshold of EW> lA, 
the fraction of galaxies with Ha emission is 5% (4% if we remove galaxies with AGN activity) 
for E galaxies, 15% (12%) for SO galaxies, 54% (44%) for SO/a-Sb and 86% (80%) for Sbc-Sd 
galaxies. We have seen that SO/a-Sb and Sbc-Sd contribute nearly the same amount to the 
luminosity density. This is due to a larger number of SO/a-Sb type galaxies than that of 
Sbc-Sd (by a factor of 1.5), whereas the fraction of emitters in the former class is less by 
1.6 times than the latter; these changes in the two factors cancel and leave the luminosity 
density unchanged. In our sample, the fraction of Ha emitters is 100% for Im galaxies. A 
high fraction 86%) is also recorded for galaxies of unclassified types. Conversely, 89% of 
emitters are spiral galaxies, and they contribute to 83% of the Ha luminosity density. Note 
that elliptical galaxies are not necessarily inactive (Fukugita et al. 2003). 

Galaxies in each morphological class contain strong Ha emitters. If one makes a cut at 
high EW threshold, say at lOA, the average of EW for those emitters look nearly indepen- 
dent of morphology. The most important difference that varies with morphology is in the 
increasing number of non-emission or weak emission galaxies towards early types. Therefore 
the average EW decreases to early types when the average is taken over all galaxies including 
non-emission galaxies. The 'average' emission equivalent width depends on how to set the 
threshold for emission galaxies. 

Kennicutt and Kent (1983) showed a clear correlation of EW(Ha-|-[NII]) with morphol- 
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ogy, notably the absence of strong emitters in earlier types (Sa or earlier). We consider this 
absence of emitters being due to their small sample, and in particular to the fact that bona 
fide early-type galaxies are selected before the observation. 

4.2. Ha emissivity and galcixy colours 

The natural assumption is that Ha emission galaxies, being hosts to vigorous star for- 
mation, should have bluer colours than galaxies that are bereft of Ha emission. Tresse & 
Maddox (1998) claimed, however, that there is no correlation between Ha luminosity L(Ha) 
and {V — I) colour. We also find in our sample that the correlation between Ha luminosity 
and g* — r* colour is weak (Figure 11a). Also, we do not see correlations with other colour 
indices, e.g., u* — g*. The correlation we should actually expect, however, is not between 
total Ha luminosity and colour, but between Ha emissivity per unit mass (or per broad band 
luminosity) and the colour. To see this point we plot in Figure lib the correlation of the Ha 
luminosity divided by i' band luminosity with gp — rp colour. The expected correlation is 
observed: there is a 2 mag change in Ha luminosity normalised by i' band luminosity across 
0.5 mag change in g* — r* colour, though the scatter is significant, especially for strong Ha 
emitters. 



5. Implications of our results 

We have derived the Ha luminosity function for the nearby universe and have made 
a break down into morphological types of galaxies. Our sample contain 665 emission line 
galaxies, which are compared, for example, to 176 galaxies in Gallego ct al. (1995), 159 
galaxies in Sullivan et al. (2000), and 110 galaxies in Tresse & Maddox (1998); see Table 4, 
which gives a comparison of the present result with those of the earlier work made for low 
redshifts. Our observed range of L(Ha) is 0.4 dex deeper than Gallego et al., although it is 
shallower by 0.4-0.8 dex than surveys of distant galaxies by Sullivan et al. (2000) and Tresse 
& Maddox (1998). 

QQ Q1 _|_r\ f\A ~l~0.10 -J Q 

The local Ha luminosity density we have obtained is 10 ' ' -0 07 h erg s" Mpc~ . 
The central value is lower than that of Gallego et al. (1995) but only by 0.08 dex, although 
this agreement is fortuitous. 

The significant discrepancy between Tresse & Maddox (1998) and Sullivan et al. (2000) 
for analyses at non-zero redshift hinders us from drawing a definitive conclusions as to the 
evolution of Ha emissivity. In comparing the characteristic luminosities we note that L* 
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and a are strongly correlated. Assuming an empirical law, dlogL* /da ^ —0.4, Tresse & 
Maddox's value is (41.64 ± 0.13 - 2 log h) dex and Sullivan et al.'s is (42.03 ±0.14-2 log h) 
dex if wc adjust a to our —1.43. Our value is in between the two. Our faint end slope is 
also between the two groups of authors. On the other hand, the Ha luminosity density of 
the two groups (io39-66±o.o4^ -|^g39.43±o.06-)^ gj,g g-i]y[pc"'^ are larger than our value by 2.2—1.3 
times. We emphasize, however, that the samples at non-zero redshifts are too small, and 
the uncertainty in the correction of the spectroscopically obtained flux to total flux might 
also be significant; the definitive conclusion may not be drawn as to the evolution at low 
redshift. The evolution can be concluded, however, if we compare our results with those at 
high redshifts (e.g., Glazebrook et al. 1999; Tresse et al. 2002). 

The conversion of Ha luminosity density to the global star formation rate per unit 
volume depends further on models of star formation and absorption of the ionising flux in 
the star-forming regions (Chariot & Longhetti 2001). If we adopt the conversion factor of 
Glazebrook et al. (1999) for the case BC96(kl96) Z/Zq = 1.0, and Salpeter IMF (with the 
lower mass cutoff at O.IM©): L(Ha) = 1.35 x lO^^erg s~^ for star formation rate of 1 Mq 
yr~^, we obtain the global star formation rate ip: 

t/j ~ O.O15/^M0yr-^Mpc"^ (7) 

ignoring internal absorption of the ionising flux. 

Some early type galaxies, even elliptical galaxies, emit strong Ha by star formation 
activity, almost as strong as that of late spiral galaxies. This contradicts the conventional 
wisdom that all ellipticals formed in the early universe. The fraction of Ha emitters, however, 
is quite small: it is only 1% if we take a threshold of EW(Ha) >10A (and 5% with EW(Ha) > 
lA): the majority of elliptical galaxies no longer have significant star formation. On the other 
hand, there are many late spiral galaxies that do not show Ha emission (in our sample all Im 
galaxies are Ha emitters). The presence of non star-forming spiral galaxies is consistent with 
the view that star formation in disc galaxies is intermittent, but the fraction of star forming 
galaxies with > IMc^yr"^ (?» 50%) indicates that in spiral galaxies such star formation 
activity takes place about half the time. The quantity that shows a strong correlation with 
morphology of galaxies is the fraction of Ha emitters, or the mean Ha luminosities, not the 
Ha luminosities of individual galaxies. 
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Table 1: Nuiiil)er of galaxies in our saiiii)lo 





T=-l ~ 6 


T=0 ~ 1 


T=1.5 ~ 3 


T=3.5 ~ 5 


T=5.5 ~ 6 


Total 


1482 


597 


518 


350 


10 


Ha detected 


665 


67 


282 


300 


10 


AGNs by Eq.(3) 


89 


13 


54 


21 






Table 2: Luminosity Function Parameters 





logL* 






log>C(HQ;) (h erg 


; s"^ Mpc"^) 


Morphology 


{h^'^ erg s^^) 


a 


(O.OOl/i^Mpc-^) 


Integrated 


Sum 


-1 < r < 6 


41.68 ±0.10 


-1.43 ±0.10 


2.78 ± 1.13 


39.31 ± 0.04 


39.33 


< r < 1 


41.02 ±0.14 


+0.79 ± 0.77 


0.78 ± 0.34 


38.13 ±0.22 


38.24 


1.5 < T < 3 


41.70 ±0.13 


-1.40 ±0.15 


1.03 ±0.43 


38.88 ± 0.03 


38.91 


3.5 < T < 5 


41.71 ±0.21 


-1.53 ±0.21 


0.96 ±0.50 


38.96 ±0.08 


39.01 


5.5 < T < 6 


- 42.76 


~ -1.77 


~ 0.006 


~ 38.13 


37.98 


T = -1 










37.80 
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Table 3: Summary of Systematic Errors 



Item 




/~^^ • /» /tt \ 

Change m L{tia) 
(factor) 


EW(Hq;) threshold (l.OA) 






decreased to —0.8 A 




1.04 


increased to 7A 




0.90 


Stellar absorption correction ±0.5A 




0.97—1.07 


entirely neglected 




1 Ad 


Extinction law (Seaton 1979) 






O'Donnell (1994), = 3.1 ± 0.2 


* 


0.98-1.10 


Errors in the estimate of Ay = 1.03 


* 


0.90-1.10 


Remove weak H/? 'AGN's 


* 


0.95 


Include all AGNs 




<1.18 


Nebulae Temperature 6000-15000 K 


* 


0.97-1.03 


Aperture corrections: 






RD94 profile with ±la 


* 


0.91-1.19 


scale with the broad band flux 




0.75-0.80 


Total systematics 




0.85-1.25 



*) Counted as systematic errors, and added in quadrature to estimate the total systematic 
errors. 



Table 4: Comparison with the earher work {h = 1) 





Gallego 95 


Tresse-Maddox 98 


Sullivan 00 


This Work 


Survey area 


471.4 deg2 


500 arcmin^ 


~ 10 deg2 


229.7 deg2 


Mean redshift 


~ 0.025 


0.21 


0.15 


0.054 


Size of the sample 


176 


110 


159 


665 


logL* 


41.56 ±0.08 


41.61 ±0.13 


42.11 ±0.14 


41.68 ±0.10 


a 


-1.3 ±0.2 


-1.35 ±0.06 


-1.62 ±0.10 


-1.43 ±0.10 


log0* 


-2.3 ±0.2 


-2.09 ±0.09 


-3.04 ±0.20 


-2.56 ±0.30 


log>C(Ha) 


39.38 ± 0.04 


39.66 ±0.04 


39.43 ± 0.06 


39.31 ± 0.04 
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Fig. 1. — Distribution of rest EW(Hq;) for each morphological type of galaxies with the 
rest EW(Hq;)> lA and the detection at > 2.5a. The lowest panel shows the total sample. 
Those plotted near the right margin are numbers of galaxies with rest EW(Hq;)> lOOA. The 
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Fig. 2.— Relation between the rest EW(Ha) and rest EW(H/3)x[Fc(H/3)/Fc(Ha)], where 
Fc is the flux density of the continuum. The region of small EW is expanded in the lower 
panel. The lines show the flt of eq. (1). 
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Fig. 3. — H7 equivalent widths (crosses), and stellar absorption equivalent widths of H7 
(triangles) plotted as a function of Ha emission equivalent widths. 



-24- 



O 
CO 



n I I I I I I I I I I I I I I I I 

331 galaxies. 



^ o 

s 



o 

CN2 



J I I I I : I I I I I L 







n I I I I \ L 



4 



K (mag) 



Fig. 4. — Distribution of extinction estimated from Balmer line ratios for 331 galaxies with 
strong emission (see text). The line shows the mean. The arrow indicates the value derived 
from eq. (1). 
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Fig. 5. — Relation between the Ha flux and the i' band flux. The dashed hues show galaxies 
measured by RD94 at varying apertures (the point moves left/upwards as the aperture 
size increases). The circles indicate the 'total flux' given by RD94. The crosses are SDSS 
spectrophotometric data for the 3" aperture, calibrated by photometrically measured fluxes. 
The dotted lines show 6 profiles that are outliers. 
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Fig. 6. — Growth curve of the Ha flux against that of the / band flux. The thick hne is the 
composite growth curve, and error bars are one sigma of the mean at several representative 
points. The thin dashed curves are those of individual galaxies in RD94. 
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Fig. 7. — Ha luminosity as a function of redshift after the aperture correction. 
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Fig. 8. — Ha LF for each morphological type of galaxy. The data points are step wise LF, 
and the curves show Schechter function fits. The histograms are actual number of galaxies 
used to derive the LF. Data symbols are: total (cross), E-SO (open circle), SO/a-Sb (solid 
triangle), Sbc-Sd (open square), Im (sohd circle). 
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Fig. 9. — Schcchter function parameters of Ha LF: (a) total sample, (b) E-SO, (c) SO/a-Sb, 
and (d) Sbc-Sd. The two contours show 1 a and 2 a errors. The cross in each panel shows 
the position of the parameter for the total sample. 
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Fig. 10. — Comparison of the rest EW distribution of Gallego et al. (1995) 's sample (dotted 
lines) with ours (solid lines). 
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Fig. 11. — (a) Ha luminosity (in magnitude scale) plotted as a function of g* — r* colour 
(in Petrosian magnitudes), (b) Ha luminosity normalised by i' band luminosity plotted as 
a function of g* — r* colour. The squares are the mean of the data and the error bars stand 



